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Abstract 

An exact numerical calculation of neutrino emissivity of two and three flavour quark 
matter have been carried out. We find that the neutrino emissivity obtained from the 
Iwamoto formula is in qualitative agreement with our calculation for two flavour quark 
matter. For three flavour case, on the other hand we find that the Iwamoto formula over- 
estimates the numerical values by ~2 orders of magnitude or more for d decay and agrees 
with the s decay results within a factor of 3-4. The dependence of the emissivity on temper- 
ature, strong coupling constant and baryon density is also quite different from the Iwamoto 
formula. 
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It has been conjectured that dense stars may consist of quark matter or quark matter 
core with neutron matter outside ( Bahcall & Wolf 1965, Burrow 1980, Freedman & McLerran 
1978, Maxwell et al. 1977, Witten 1984 ). Although theoretical understanding of the properties of 
quark matter is not yet available, various quark models have been used to calculate the equation 
of state of the quark matter and determine the properties of quark stars ( Farhi & Jaffc 1984, 
Ghosh & Sahu 1993, Goyal & Anand 1990, Haensel et al. 1986 ). Unfortunately, it is found 
that the properties of quark stars, such as surface gravitational redshift z,moment of inertia I, 
maximum mass M, radius R and pulsar periods P, are not significantly different compared to 
those of neutron stars. Therefore it is difficult to distinguish one from the other observationally. 

On the other hand Iwamoto ( Iwamoto 1980, 1982 ) has proposed that neutrino emissivity 
( e ) could play a significant role in distinguishing between quark and neutron stars because 
it differs by orders of magnitude for the two. Particularly e for quark stars is larger by 6-7 
orders of magnitude than neutron stars which could lead to faster cooling rate for quark stars, 
thus reducing their surface temperature. There are, however, a number of other mechanisms ( 
pion condensation ( Brown 1977 ) and modified URCA processes ( Friman & Maxwell 1979 ) ) 
proposed to increase e for neutron matter. 

Iwamoto ( 1980 ) has derived the formula for e using apparently reasonable approximations 
and this formula has been widely used ( Alcock et al. 1986, Datta et al. 1988, Duncan et al. 
1983, Goyal & Anand 1990 ) to calculate e for two and three flavour quark matter. According 
to his formula e is proportional to baryon density { Ub ) , strong coupling constant ( etc ) and 
sixth power of temperature (T) for d quark decay. For s quark decay T dependence of e is same 
as that for d decay. Furthermore his results show that electron and quark masses have negligible 
contribution to e and s quark decay ( in case of three flavour quark matter ) plays a rather 
insigniflcant role. 

In the present letter we want to report an exact numerical calculation of e and a comparison 
of our results with the Iwamoto formula. Our results show that the Iwamoto formula is in 
qualitative agreement with our calculation for two flavour quark matter, although it yields e 
which is 3-4 times larger. However we flnd that the Iwamoto formula overestimates e by about 
2 orders of magnitude or more for three flavour quark matter and hence it should not be used 
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for three flavour quark matter. In the following, we first discuss the equilibrium composition of 
matter then we sketch the calculation of e and present the results. 

The simplest possible processes for neutrino emission in two flavour degenerate quark matter 
are direct (5 decay reactions 

d — > + e~ + Pg 

u + e~ ^ d + Ve (1) 
The chemical equilibrium and charge neutrality implies 

and 

2nu -Ud-Sue^O (3) 

The baryon density is defined as ub = {nu + nd)/3 where rii = g-P%J {Qt^"^) is the particle number 
density. Degeneracy factor g is 6 for quarks and 2 for electron. 

For three flavour degenerate quark matter the neutrino emission occurs via d as well as s 
decay. The s decay reactions are 

s — >• w + e~ + Pe 

u + e~ ^ s + (4) 
Prom equations (1), (2) and (4), the chemical equilibrium and charge neutrality gives 

IJ'd = l^s (5) 

and 

- Ud - Us - Sue = (6) 

Baryon density is defined as ub — {riu + nd + ns)/S. Thus, given the baryon density, quark 
and electron number densities and chemical potentials are calculated using the constraints above. 
These are then used to evaluate neutrino emissivity. 



Neutrino emissivity ed(s) for d (s) decay is given by ( Iwamoto 1980,1982 ) 

^d(s) - ^d{s) / " Pd{s)<i Pud Ped Pu — ^tt; ^ 

x^'^iPdis) - Pu - Pe - Pv)n{pd{s))[^ - n{pu)\[l - n{pe)\ (7) 

where Pd{s)-i Pm Pe, Pv are the four momenta of d(or s), u, e, v] n{pi) {i— u,d (s), e) is the 
Fermi distribution given by 

and 

where G is Weak coupling constant and Qc is Cabibbo angle, Ei is the energy and is the 
momentum. In terms of Fermi momentum pp the quark chemical potential \x ( upto lowest order 
in etc ) is ( Baym & Chin 1976 ) 

= - + 1 ln{x + 77 11 

X 37r xri 

where x = pF/m and r] = ^J(1 + x^), m being the quark mass. For massless quarks eq(ll) 
reduces to 

= (1 + (12) 

Since in our calculation we have taken u, d quark masses to be zero and s quark mass to be 
nonzero ( 150 - 200 MeV), the energy momentum relation of quarks is approximated by eq(ll) 
for s quarks and eq(12) for u and d quarks, that is, ^ is replaced by E and Pf is replaced by p 
in above eq(ll-12). This is reasonable as only the energies and momenta close to Fermi surface 
contribute to the matrix element. For massless electrons E^, = p^ and the energy delta function 
is used to perform pe integral. This gives 

^ {Ed(s) - Euf - (Pdjs) - Pnf - 2pd(s)P„(l - COSe^) 

{2p^cose^e + '2pd{s)Cos9, + 2Ed^s)-'2E^) ^ ' 
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where 6u, Og and Que are the angles between d and u, d and e and u and e respectively. The 
momentum delta function in eq(7) is used to integrate over the neutrino momentum p,^. Of the 
rest of the integrals dQ.dd(t>u integral gives Stt^ and we get 



In the above expression (eq(14)) limits for momentum integrals are from to oo, for cosOu, 
cosOe from -1 to +1 and for 0e from to 27r . However we find that the integrand is a sharply 
peaked function about p„ ^ Pf{u), pa ~ Pf{d) (or ~ P/(s)) , 0e ~ and certain values of 
cosOe and cos9u{< !)■ Thus one must ensure that there are sufficient number of integration 
points in this region. This we have done. We have also ensured the convergence of the integral 
with respect to the number of integration points and step length for integration. 

Two and three flavour e, for different T, ac and hb arc given in Tables 1 and 2 (a, b) 
respectively. Corresponding Fermi momenta are given in Tables 3 and 4. For two flavour quark 
matter varies as which is slightly lower than the Iwamoto results . In fact the T 

variation changes with temperature as well as etc and ub- We also find that dependence of e 
on etc and ub are approximately e oc Ub^^ and e oc a'^^^ respectively. These are close to values 
obtained by Iwamoto and others ( Iwamoto 1982, Alcock et al. 1986, Datta et al. 1988, Duncan 
et al. 1983 ) . We find that with increase in baryon density, electron fraction in two flavour matter 
increases, which seems to play a crucial role resulting in the increase in e. The calculated e from 
the Iwamoto formula has also been displayed in Table 1. A comparison shows that Iwamoto 
results are larger by factor of 3 or more with the difference increasing with T and ub- Thus we 
find that Iwamoto formula provides a reasonable approximation to an exact calculation if the 
difference of a factor of 3-4 is acceptable. 

The neutrino emissivity and its dependence on T, etc and ub are entirely different for three 
flavour quark matter. For d decay varies as T'3-6-4.7 f^j. _ i^QMeV and T^-^^ for = 
200MeV. The etc and ub dependences also vary in a wide range; for T=0.1 MeV oc and 
Ed (X n^^, where as for T= 0.6 MeV oc and oc n^^, being 150 MeV. Unhke the 




(14) 
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two flavour case, here electron density decreases with increasing baryon density. Corresponding 
€4 also decreases. The ea obtained using the Iwamoto formula comes out to be much larger than 
our values as shown in Table 2a and 2b. For s decay e has been found to be much different when 
compared with d decay. We find e, oc T^ ^-s-ss foj. ^ 150MeV and T^-^-^-^ for m, = 200MeV. 
Further increase in results in increase in power of T beyond 6. The dependence of is 
similar to that of e^. One particular change to be noticed is that for ac = 0.05, eg decreases with 
density, whereas for etc = 0.1 it increases with density. This difference in s decay compared to d 
decay implies the significant role played by rris. 

Here we would like to note that the cases where Iwamoto formula agrees reasonably well 
with our result, the difference between Pf{u) + pf{e) and Pf{d) ( or Pf{s)) is much larger than 
the temperature. On the other hand, when this difference is smaller or comparable with the 
temperature, the Iwamoto formula overestimates the exact result by orders of magnitude. This 
behaviour is also seen when the temperature is larger. 

In order to investigate this further, we have attempted to evaluate the integral approxi- 
mately. For this we note that the integrand peaks at certain values of the integral variables ( 
Pd{s){'^) ■iPu{'^) 1 cos9u{m) ,cos9e{'m) and 4>e{m)). We find that a reasonable Gaussian approxima- 
tion for the integration is 



The above integration eq(15) can be performed analytically if the integration limits are from 
—00 to -|-oo and the result is 

7r5/2 

e,(,) = C.8n^ (16) 

We find that for both two and three flavour quark matter, the emissivity from above approx- 
imation (eq(16)) are within factor of 2 of our exact numerical results. Further we find that for 
two flavour case C oc T and oc ^, so that e oc (approximately). For three flavour d decay. 



5 



C is not linear in T (~ T^/^) and the dependence of a^'s on T is also different. Thus, in this 
case the temperature dependence obtained after numerical integration differs significantly from 
Iwamoto results. 

It was first pointed out by Duncan et al. ( 1983 ) that treating electron fraction as a constant 
( Yf. = 0.01 by Iwamoto ) is not correct and a detailed treatment of the equilibrium composition 
of quark matter is needed. For a fixed density they solve the chemical equilibrium and number 
density conditions to get the prii) and We indeed follow the same method. Thus the 

neutrino emissivity denoted by edi{si) in Table 2 ( a, b ) are essentially the three flavour e 
calculated by Duncan et al. ( 1983 ). 

Our calculation shows that for 2-fiavour quark matter, e calculated from the Iwamoto formula 
agrees qualitatively with the exact calculation and therefore earlier conclusions ( Alcock et al. 
1986, Datta et al. 1988, Duncan et al. 1983, Iwamoto 1982 ) regarding cooling rates of 2-flavour 
quark matter are still valid. However for 3-flavour quark matter e is smaller than the results 
obtained using the Iwamoto formula by 2-3 orders of magnitude or more for d decay and factor of 
3 or more for s decay. As a result, e of three flavour quark matter is atleast 2 orders of magnitude 
smaller than that predicted by the Iwamoto formula. Therefore, the cooling rates of three flavour 
quark matter are considerably smaller than the rates estimated from the Iwamoto formula. 

Our calculation shows that e of two flavour quark stars is about 3 or more orders of magnitude 
larger than that of three flavour quark stars. Thus, a two flavour quark star will cool more rapidly 
than a three flavour quark star. On the other hand the dynamics of the quark matter suggests 
that the three flavour quark matter is most stable ( Witten 1984 ) and a two flavour quark matter 
will decay into three flavour quark matter by weak interactions. We are not aware of calculations 
where this decay rate has been calculated. But, as the two to three flavour conversion rate as well 
as neutrino cooling involve weak interactions, the two rates would be comparable. Furthermore, 
there are arguments ( Alcock et al. 1986 ) to suggest that quark stars, when formed, may 
consist of three flavour quark matter. Therefore neutrino cooling rates of quark stars may be 
much smaller than those suggested from earlier calculations. In this context, we feel that the 
knowledge of two to three flavour decay rate will be useful. 
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Table 1. Neutrino emissivity for two flavour quark matter. Here e^/ is the emissivity calculated 
using Iwamoto formula. 



T 


riB 




-- 0.1 




0.05 


(MeV) 








Cd 








{erg/cm? / s) 


[ergjcw? 1 s) 


{erg/cm? / s) 


[ergjcnv' j s) 




0.6 


9.93x10+25 


3.01x10+26 


5.21x10+25 


1.45x10+26 


0.1 


1.0 


1.68x10+2^ 


5.01x10+26 


8.86x10+25 


2.41x10+26 




1.4 


2.35x10+26 


7.02x10+26 


1.25x10+26 


3.38x10+26 




0.6 


3.35x10+29 


1.23x10+3° 


1.15x10+29 


5.93x10+29 


0.4 


1.0 


5.84x10+29 


2.05x10+30 


2.28x10+29 


9.88x10+29 




1.4 


8.34x10+29 


2.87x10+30 


3.50x10+29 


1.38x10+30 




0.6 


3.12x10+=^° 


1.40x10+31 


7.50x10+29 


6.75x10+30 


0.6 


1.0 


5.77x10+^° 


2.34x10+31 


1.62x10+30 


1.12x10+31 




1.4 


8.50x10+3° 


3.27x10+31 


2.65x10+30 


1.57x10+31 



8 



Table 2a. Neutrino emissivity for three flavour matter. Here e^/ and Sgi are the emissivity 
calculated using Iwamoto formula. 



ac 




T 


riB 












(MeV) 


(MeV) 




((7Y/ / Clll^ /s) 


(cr(j/ciii-'/.s) 


{cr(j/ ciii^ /s) 


(crg/cin^/s) 








0.6 


1.31x10+23 


8.08x10+24 


7.47x10+23 


1.73x10+24 






0.1 


1.0 


2.11x10+22 


5.73x10+24 


4.81x10+23 


1.23x10+24 








1.4 


3.46x10+21 


3.56x10+24 


2.34x10+23 


7.60x10+23 








0.6 


3.19x10+25 


2.31x10+28 


2.23x10+27 


7.10x10+27 




150.0 


0.4 


1.0 


9.73x10+24 


2.35x10+28 


8.15x10+26 


5.03x10+27 








1.4 


4.21x10+24 


1.46x10+28 


1.62x10+26 


3.11x10+27 








0.6 


2.07x10+26 


3.77x10+29 


1.84x10+28 


8.09x10+28 






0.6 


1.0 


8.48x10+25 


2.67x10+29 


5.02x10+27 


5.74x10+28 


0.1 






1.4 


4.75x10+25 


1.66x10+29 


8.74x10+26 


3.55x10+28 








0.6 


2.87x10+24 


2.39x10+25 


2.39x10+24 


5.13x10+24 






0.1 


1.0 


1.31x10+24 


2.17x10+25 


2.18x10+24 


4.66x10+24 








1.4 


6.02x10+23 


1.93x10+25 


1.86x10+24 


4.14x10+24 








0.6 


7.98x10+26 


9.79x10+28 


8.42x10+27 


2.10x10+28 




200.0 


0.4 


1.0 


2.88x10+26 


8.90x10+28 


7.41x10+27 


1.91x10+28 








1.4 


1.30x10+26 


7.90x10+28 


6.09x10+27 


1.69x10+28 








0.6 


3.77x10+27 


1.12x10+^° 


9.30x10+28 


2.39x10+29 






0.6 


1.0 


1.47x10+27 


1.01x10+30 


7.71x10+28 


2.17x10+29 








1.4 


7.44x10+26 


9.00x10+29 


5.70x10+28 


1.93x10+29 
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Table 2b. Neutrino emissivity for three flavour matter. Here eai and egi are the emissivity 
calculated using Iwamoto formula. 





rus 


T 














(MeV) 


(MeV) 


{.for') 


((7-.(//r7;/'V'''-) 


{crcj/ciii^'/.s) 


(crg/ciii^'/s) 


(crg/ciii^/s) 








0.6 


4.34x10+23 


1.14x10+25 


1.69x10+24 


0.00x10+00 






0.1 


1.0 


2.84x10+23 


1.26x10+25 


2.18x10+24 


5.29x10+24 








1.4 


2.09x10+23 


1.34x10+25 


2.40x10+24 


5.62x10+24 








0.6 


6.73x10+25 


4.68x10+28 


7.79x10+27 


0.00x10+00 




150.0 


0.4 


1.0 


4.26x10+25 


5.16x10+28 


8.75x10+27 


2.17x10+28 








1.4 


3.15x10+25 


5.48x10+28 


9.21x10+27 


2.30x10+28 








0.6 


3.01x10+26 


5.33x10+29 


8.79x10+28 


0.00x10+00 






0.6 


1.0 


2.01x10+26 


5.88x10+29 


9.46x10+28 


2.47x10+29 


0.05 






1.4 


1.56x10+26 


6.24x10+29 


9.64x10+28 


2.62x10+29 








0.6 


2.90x10+24 


2.27x10+25 


3.08x10+23 


0.00x10+00 






0.1 


1.0 


2.27x10+24 


2.53x10+25 


2.29x10+24 


0.00x10+00 








1.4 


1.85x10+24 


2.70x10+25 


3.76x10+24 


0.00x10+00 








0.6 


6.70x10+26 


9.29x10+28 


1.09x10+28 


0.00x10+00 




200.0 


0.4 


1.0 


4.37x10+26 


1.03x10+29 


1.58x10+28 


0.00x10+00 








1.4 


3.24x10+26 


1.11x10+29 


1.81x10+28 


0.00x10+00 








0.6 


2.84x10+2^ 


1.06x10+30 


1.47x10+29 


0.00x10+00 






0.6 


1.0 


1.86x10+27 


1.18x10+30 


1.89x10+29 


0.00x10+00 








1.4 


1.39x10+27 


1.26x10+30 


2.10x10+29 


0.00x10+00 
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Table 3. Baryon number density Ub, Fermi momenta of u-quark Pf{u), d-quark pr^d) and 
electron ^^(e) for different where — Pf{u) + Pf{^) — PF^d) . 





riB 


Pf{u) 


PF{d) 


PF{e) 


Apd 






(McV) 


(MeV) 


(MeV) 


(MeV) 




0.60 


357.80 


449.20 


99.15 


7.75 


0.1 


1.00 


424.22 


532.58 


117.56 


9.20 




1.40 


474.57 


595.79 


131.51 


10.29 




0.60 


357.71 


449.25 


95.43 


3.89 


0.05 


1.00 


424.11 


532.65 


113.15 


4.61 




1.40 


474.45 


595.87 


126.57 


5.15 
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Table 4. Baryon number density n^, Fermi momenta of u-quark Pf{u), d-quark prid), 
Pf{s) and electron pF{e) for different and different etc , where Ap^ — Pf{u) +Pf{g-) 
and Aps = Pf{u) +PF(e) - Pf{s) 



s-quark 
-pF{d) 



rris 


etc 




Pf{u) 


PF{d) 


Pf{s) 




Apd 


Aps 


(MeV) 






(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 






0.60 


356.99 


360.06 


353.86 


3.33 


0.26 


6.46 




0.1 


1.00 


423.26 


424.81 


421.69 


1.69 


0.14 


3.26 


150.0 




1.40 


473.49 


474.27 


472.72 


0.84 


0.06 


1.61 






0.60 


356.99 


365.89 


347.62 


9.28 


0.38 


18.65 




0.05 


1.00 


423.26 


430.29 


415.98 


7.33 


0.30 


14.61 






1.40 


473.49 


479.49 


467.34 


6.25 


0.25 


12.40 






0.60 


356.99 


365.93 


347.58 


9.70 


0.76 


19.11 


200.0 


0.1 


1.00 


423.26 


429.10 


417.25 


6.34 


0.50 


12.35 






1.40 


473.49 


477.66 


469.25 


4.52 


0.35 


8.76 






0.60 


356.99 


374.26 


337.85 


18.00 


0.73 


37.14 




0.05 


1.00 


423.26 


437.14 


408.39 


UA7 


0.59 


29.34 






1.40 


473.49 


485.45 


460.90 


12.46 


0.50 


25.05 
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